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Abstract 

The static as well as the dynamic behaviour of granular material are de- 
termined by dynamic and static friction. There are well known methods to 
include static friction in molecular dynamics simulations using scarcely un- 
derstood forces. We propose an Ansatz based on the geometrical shape of 
nonspherical particles which does not involve an explicit expression for static 
friction. It is shown that the simulations based on this model are close to 
experimental results. 
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The behaviour of fluidized dry granular material, like sand or powder, reveals a rich 
variety of effects, which can not be observed in other substances. Those effects have been 
observed and investigated by many scientists over a long period [[|. Examples for the 
most interesting effects are fluidization, convection cells and heap formation under vibration 
@F U> s i ze segregation ("Brazil-nut" effect) [JTtJ- [0 deformation under shear force |I3 



shape segregation of differently shaped grains in a pipe |L4| and clustering instabilities [|15 



Density waves emitted from outlets |Tg inside material flowing through pipes [17|] and at 



the surface of an inclined chute |L8[ have been intensively investigated. Recent results gave 
evidence that convection cells due to walls or amplitude modulations play an important role 
in the process of the formation of macroscopic structures 0. Of particular interest are the 
dynamic as well as the static behaviour of avalanches going down the slope of a sand pile. 
Theoretical as well as experimental investigations |]T9]] — [^3| led to the hypothesis that their 



mass and their time distributions can be described by the self organized criticality-model. 
There are experiments, however, that do not agree with this hypothesis || [p^j . Recently 
many experimental observations have been reproduced by numerical simulations. There is a 
wide variety of simulation methods including Monte-Carlo simulations (e.g. [fCT]]), molecular 
dynamics simulations (e.g. || |TD[ |T^]), and random walk approaches [F25J. These simulations 



gave many interesting information on the microscopic effects underlying the behaviour of 
macroscopic amounts of granular material. Many of the effects observed in experiments are 
consequences of static friction between the grains. In most of the current simulations special 
terms for static friction are used to mimic static behaviour of granular material e.g. j|] p6| . 
The aim of this paper is to show that it is possible to reproduce the experimental results 
by molecular dynamics simulations without introducing such a static friction force but by 
simulating nonspherical particles. We show that our simulations with nonspherical particles 
agree better with experimental results than equivalent simulations introducing static friction 
forces as it is usually done. 

Since it is extremely complicate to calculate collisions of cubic particles we choose in two 
dimensions particles similar to squares but consisting of spheres. A further advantage of 



this model is that we are able to vary the shape steadily from a sphere almost to a square. 
A related Ansatz for nonspherical grains was recently done by Gallas and Sokolowski |27| . 



there each grain consists of two spheres rigidly glued to each other. Each of our nonspherical 

(k) 

particles k consists of four spheres with equal radii r\ , located at the edges of a square 
of size Lft\ and one sphere with radius r!$ = — rf^ in the middle of the square 

(fig. [I]). For the case that two spheres % and j of the same particle k or of different particles 
touch each other during a collision there acts the force 



c 



F 



Y(ri + Tj - \xi - Xj\) + rym e ff\xi - Xj\ 



. rrii ■ rrij 

With m eff ~ 



rrii + rrij 

where Y is the Young modulus and 7 is the phenomenological friction coefficient. In addition 
to the forces between each two particles of the system, there are forces between each pair of 
spheres i,j where i and j both belong to the same grain, due to a damped spring 



F s 



a(C {k) - \xi -xj\) + ^sp^lxi ~ ■>-, 



where a and 7s p are the spring constant and the damping coefficient. If the spheres % and 
j are both located at the same edge of the square then equals lS k \ if i is the central 
sphere then = L^/y/2. 

The dynamics of large numbers of such nonspherical particles was investigated simulating 
the flow of granular material in a rotating cylinder under gravity. For the integration we 
used a sixth order predictor-corrector method [p8|| . In a cylinder of diameter D = 260 



we simulated the flow of 1000 nonspherical particles of different size L^ k > with Gaussian 
probability distribution p(L^) = N(3, 1) each consisting of five spheres. For the parameters 
we chose Y = 10 4 kg/s 2 , 7 = 1.5-10 4 s~\ a = 10 4 kg/s 2 , 7 5p = 3-10 4 s^andrf = 1/4-L^. 
The cylinder consists of spheres with different radii to mimic a rough surface. The mean 
value of these spheres equals the mean value of the L^ k >. The cylinder was rotated clockwise 
with very low angular velocity Q. During the uniform rotation of the cylinder the flow of 



the grains was very inhomogeneous, due to avalanches going down the inclined surface. This 
behaviour is called stick-slip motion. The time evolution of the slope of the surface as 
well as the averaged velocity v of the particles at the surface for a fixed angular velocity 
Q = 0.002 s -1 are drawn in fig. |2] (curves v (a), 9(a)). The angle was plotted in rad, v 
in 50 • Since the number of particles is not too large our surfaces are not smooth. 

Hence we have to determine the inclination indirectly as the angle between the straight 
line connecting the centre of mass point of the grains and the middle point of the rotating 
cylinder and the direction of gravity. The angle and particularly the average velocity of the 
surface particles fluctuate drastically and irregularly as it is typical for stick-slip motion. 
This behaviour was observed experimentally before by Briscoe, Pope and Adams J29| and 
by Rajchenbach []30[| . The plots v(b),9(b) in the same figure show the equivalent data for 
the simulation using spherical particles. The radii of the spheres were Gauss-distributed too 
with p(ri) = N[l, 1]. The spherical grains undergo the same force as the spheres of which the 
nonspherical particles consist. To mimic static friction we include for the case of spherical 
particles rotation as a further degree of freedom of the grains and add the force 

T*^f — vyi in J r\/ <m „ „ I nt . I 1 1 \ TT 1 ^ I 



F-f = min{-7 s m e// |tf reZ |,/i|i^.|} ( 



J I iXj ^ iic j | 



with v re i = (xi — Xj) + Ti ■ uJi + Tj ■ ujj , 

where cJj is the angular velocity of the i-th particle, 7 S is the shear friction coefficient, and \x is 
the Coulomb parameter (j a = 3-10 4 s _1 , \i = 0.5). This Ansatz is the most popular to include 
static friction between particles which roll on each other into the expressions for the forces 
used in molecular dynamics simulations. It was introduced by Cundall and Strack and 
modified by Haff and Werner [M. Most of the molecular dynamics simulations of granular 



material base on this Ansatz. The force F*f was implemented only for the simulation of 
spherical grains but not for the nonspherical. 

Obviously the qualitative behaviour of the slope G resemble each other in both simu- 
lations but quantitatively we get for nonspherical grains more than twice the mean angle 



(O ns ) than for spherical (0 sp ). For very low rotation velocity Q = 2 ■ 10 -3 we found O sp = 7° 
and O ns = 19°. In the experiment |3(| was measured O ~ 27°. The average velocity of the 
surface grains differs significantly too for both cases. The typical avalanches in the case of 
nonspherical particles can not be observed for spheres. The curve v (b) is much smoother. 
In the experiment one observes stick-slip motion ||30|| . Fig. |3] shows the slope O of the surface 
as a function of the angular velocity of the cylinder Q for both nonspherical and spherical 
grains. In both cases the curves are close to a straight line. For much higher angular veloc- 
ities than used in our simulations the grains do not move stick-slip like but continuously. 
In this regime was experimentally found tt ~ (0 - c ) m , with m = 0.5 ± 0.1 Pfl. With 
the same Ansatz we find m ~ 1.25 for the stick-slip regime. As shown above the simulation 
with nonspherical grains coincides much better with the experimental observations than 
equivalent simulations using spheres. 

In our second simulation we investigate the evolution of a stable pile of granular material 
by continuously dropping particles on the top of the pile due to the experiment of Held 
et. al. [plj] . Beginning with an empty rough plane we drop the next particle when the 
maximum velocity vanishes, i.e. when it is smaller than a given very small threshold v max <C 
1. The rough plane was simulated by a chain of fixed spheres of random radii with mean 



T~i = LW where is the size of the z-th nonspherical grain. The parameters of the 
simulation were the same as in the previous experiment. 

During the simulation we noticed that the slope of a pile of nonspherical grains does 
not depend on the number of particles. For spherical grains, however, the heap dissolves 
under gravity, i.e. the slope decreases with increasing particle number. There are molecular 



dynamics simulations of stable piles with spherical grains, e.g. pq| , but there the particles 
are not allowed to roll on each other, hence they can only slide, this behaviour does not 
correspond to experimental reality. 

If the platform above which the heap is built up has a finite length P one can investigate 
the fluctuations of the mass rrih of a heap of definite size and the distribution of the size of 
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the avalanches, i.e. the mass fluctuations of the heap. In fig. f| is drawn the time series of 
the mass rrih for fixed P. The mass fluctuates irregularly due to avalanches of different size 
going down the surface of the heap. The size distribution of the avalanches follow a power 
law, fig. | shows the log-log-plot of the spectrum. For the exponent H{Na) ~ {Na) 171 we 
found m ~ —1.4. The experiments yield m ~ —2.5 [21] and m m —2.134 |2"2fl . For the case 
of spherical grains we cannot find avalanches. 

The ratio between the size of a grain and the radii of the spheres at the corners determines 
whether the grains shape is closer to a sphere or to a square. Hence we define a shape value 



S — 1 R^in/R 



cc 

max > 



where R^ in and R^ ax are the extremal values of the distance between the convex cover of 
the nonspherical grain and its central point (fig. |6]). For the limit S — > the grains have the 
shape of spheres. The function reaches its maximum S m = 0.255 for a grain which convex 
cover is most similar to that of a square. To investigate the influence of the the shape S 
of the grains on the result of the simulations we have to scale the density p of the material 
the grains consist of to ensure that the total mass of each grain remains constant. Fig. |5| 
shows the angle of the heap as a function of the shape S. For grains with shape S = S m , 
which corresponds to /h )s m — 9.66 , the inclination of the heap reaches a maximum 
too. The angle $ ~ 23.1° agrees with experimental data, Bretz et. al. p2 found $ ~ 25°. 



Each other value S corresponds to two different particle shapes both closer to a sphere than 
the S'm-particle. The values marked by are due to grains with L^/rf ] < (L^/r ( i k) ) Sm , 
+ designates the slope of the heap for particles with L^/r^ > (l/ fc ) /r!j®)g m . As expected 
the slope $ of the heap rises with growing L( fc )/rf } up to S reaches its maximum S = S m . 
For larger ratios L^/r^ (S < S m ) the slope $ decreases. The dashed line in fig. || displays 
the inclination $ sp we measured for a heap of spheres, which corresponds to S —>■ 0. The 
value $ sp gives a lower boundary for the slope. The observed $ values for S G (0, S m ) lie 
between and $>(S m ). The calculation for the data shown in fig. |6| are very computer 
time consuming. For this reason we are not able to present more data. 



The simulations described above demonstrate that nonspherical grains are able to de- 
scribe the static behaviour of granular materials, such as stick-slip motion in a rotated 
cylinder at very low angular velocity and the angle of repose of a sandpile. It is shown 
that equivalent simulations with spherical grains and an additional term which describes 
the static friction due to the Coulomb law could neither reproduce the experimental results 
for stick-slip motion nor for the angle of repose of a sand pile. The angle of repose reaches 
its extremal value for grains which shape is close to a square. 

Hence we conclude that our microscopic model of nonspherical grains supplies a possi- 
ble description of the static behaviour of a granular material. The results regarding non- 



sphericity agree well with those in [27 



The authors thank J. A. C. Gallas and H. J. Herrmann for stimulating and enlightening 
discussions. 
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FIGURES 



FIG. 1. Shape of a nonspherical particle. 

FIG. 2. Time evolution of the slope 6 of the surface and the averaged velocity v of the particles 
at the surface of the flow for nonspherical (a) and spherical (b) grains. Due to avalanches 
fluctuates significantly, while in case (b) occur only very small avalanches. 

FIG. 3. The inclination of the surface as a function of the angular velocity of the cylinder. 

FIG. 4. Total mass of a pile of nonspherical grains on a platform of finite length P = 820. 

FIG. 5. Size distribution of the avalanches. The line displays the function }i(Na) = (Na)~ 1A - 

FIG. 6. Slope $ of a heap over the shape value S for grains with L^> /rf^ < {L^> /rf^)s m (©) 
and L^/r\ k) > {L^ /rf ] ) Sm (+)• The dotted line leads the eye to the function $ = 130-S+const. 
The dashed line displays the inclination observed in simulation with spherical particles. 
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